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We calculate RPA-BCS based spin resonance spectra of the newly discovered iron-selenide superconductor 
by using the two orbitals tight-binding model in 1 Fe unit cell. The slightly squarish electron pocket Fermi 
surfaces at (tt, 0)/(0, tt)— momenta produce leading interpocket nesting instability at incommensurate vector 
q ~ (tt, O.Stt) in the normal state static susceptibility, pinning a strong stripe-like spin-density wave or antifer- 
romagnetic order at some critical value of U. The same nesting also induces d^2 _y2 —pairing in one Fe unit cell. 
The superconducting gap is nodeless and isotropic on the Fermi surfaces as they lie concentric to the four-fold 
symmetric point of the d— wave gap maxima, in agreement with various experiments. This produces a slightly 
incommensurate spin resonance with upward dispersion, in close agreement with neutron data on chalcogenides. 
Finally, we demonstrate the conversion procedure from a 1 Fe unit cell to a 2 Fe one in which the gap symmetry 
transformed simultaneously into a da;y— pairing and the resulting resonance spectra moves from the q ~ (tt, tt) 
to the q ~ (27r, 0)/(27r, 0) region. 

PACS numbers: 74.7().-b.74.20.Rp.74.20.Pq,74.20.-z 



I. INTRODUCTION 

The newly discovered high-Tc superconductor in iron- 
selenide (Fe2Se2-) based compounds' not only increases the 
number in the unconventional superconducting (SC) family, 
but also bears on several unique magnetic and SC properties 
coming from its unusual Fermi surface (FS) topology. Both 
angle-resolved photoemission spectroscopy (ARPES)^ and 
the local density approximation (LDA) calculations^^, have 
demonstrated that in Tlo.63Ko.37Fe2_2:Se2 at a critical value 
of Fe vacancy x — 0.22 to tune the superconductivity, one 
can eliminate the hole-pocket at the F point and the FS con- 
sists of only two electron pockets at (tt, 0)/(0, tt). Specific 
heat'', and ARPES"*^ studies find that these FSs host node- 
less and isotropic SC gaps, whereas NMR measurements pre- 
dict that the corresponding pairing is spin-singlet.^^ Evidence 
is mounting from Raman," /xSR,'^ optical,'^ NMR,'SIII1 and 
other experiments"', and supported by LDA calculationJ^H 
that these materials possess a magnetic order ground state 
which coexists with superconductivity. 

In cuprates, the FS is a hole-like pocket centering at the 
k = (tt, tt) point and the nesting between inter FSs is strongest 
along Q ~ (tt, tt) which is responsible for a spin-density wave 
(SDW), (i J. 2 _ J, 2— wave pairing, and a spin resonance peak at 
Q with an 'hour-glass' dispersion in the SC state! 1 In a 
Ce-based heavy fermion, the hole pockets at fc = (tt, tt) lead 
to Q = (tt, tt, tt)— resonance with dj.2_y2— wave symmetry.'^ 
In iron-pnictide^ and iron chalcogenides,^ '^^"^ the FS has 
two hole-pockets at F and two electron-pockets at fc = 
(tt, 0)/(0, tt) in the unfolded 1 Fe unit cell and the leading 
nesting between the hole and the electron pockets constitute 
q — (tt, 0)— SDW, s^— pairing, and Q — (tt, tt)— spin res- 
onance with an 'hour-glass' or upward dispersion (only ob- 
served in chalcogenides). In KFe2As2 only the hole pocket is 
present at T and is predicted to have nodal d— wave pairing. 
On the contrary, in iron-selenide only one (or two concentric) 
electron-pocket(s) appeared at fc = (tt, 0)/(0, tt) at a partic- 



ular value of Fe vacancy.!^ A similar FS topology appears 
in the particular case of electron-doped cuprate near under- 
doping where the antiferromagnetic (AFM) gap eliminates the 
hole-pocket in the nodal region, resulting in a nodeless d-wave 
gap.2<' 

The role of FS topology and the interplay between electron 
and hole pockets and the exact shape of the neutron scatter- 
ing intensity as a consequence of the change in the FS shape 
across these compounds remain a major challenge in the field. 
In this paper we investigate the unusually located electron 
pockets in iron-selenide superconductors in a two band tight- 
binding (TB) model and their role in magnetic susceptibility. 
We find that the strong instabiUty in the static susceptibility 
evolves around an incommensurate vector q — (tt, O.Stt) at 
some critical value of U leading to a stripe like SDW or AFM 
order in these systems, in agreement with experiments .liiHi^ 
Unlike the stripe order in cuprates or other iron based super- 
conductors, here the stripe SDW order has the same q modu- 
lation at which a d^2_y2~ pairing symmetry possesses a sign 
change at the hot-spots which results in a spin resonance peak 
at Q ~ (tt ± (5, TT ± (5) in the SC state. Focusing mainly on 
Tlo.63Ko.37Fei.78Se2, and at the experimental SC gap value 
of 8.5 meV,^ we predict the resonance peak to be present in 
the range of 12.4meV at Q ~ 0.78(7r, tt). The resonance pro- 
file also yields an 'hour-glass' or upward dispersion and 45° 
rotation of the resonance profile, in close agreement with ob- 
servations in chalcogenide^23|22]^jj(j cuprates.G^ 

An important aspect of identifying the d— wave pairing and 
the resonance is that all these results have to be converted 
simultaneously to the actual 2 Fe unit cell in which they 
will look as follows: (1) stripe order nesting (tt, O.Stt) — >■ 
(O.Stt, O.Stt), (2) d^2_y2-wa\e pairing— > d^y-wawe and (3) 
spin resonance at (tt ± (5, tt ± (5) — > (2Tr ± 6,0). Further- 
more, all these results are solely governed by the FS topology 
and thus they are reproducible with the inclusion of a higher 
number of bands away from the FS. 

The paper is organized as follows. In Section II, we present 



our two-orbital TB formalism and the fitting of the dispersion 
and FSs to first-principle calculation and ARPES data, respec- 
tively. The evolution of the 'stripe' competing order through 
the calculation of static susceptibility within the random phase 
approximation (RPA) calculations is given in Section III. The 
spin-resonance spectra and the role of wave pairing are 
studied in Section IV. Section V gives a detailed discussion 
of how all the results of the 1 Fe unit cell Brillouin zone (BZ) 
transformed to the 2 Fe unit cell BZ. The temperature depen- 
dence of the spin-resonance spectra at Q = (tt, tt) in 1 Fe 
unit cell is studied in Section VI. Using a five band model, we 
show in Section VII that all the results obtained with two-band 
model remains unchanged as the higher orbitals lie away from 
the Fermi level. Finally, we conclude in Section VIII. 



II. TIGHT-BINDING FORMALISM 

First principle calculations show that the main contribu- 
tion to the density of states (DOSs) near Ep comes from 
i2g~orbitals of Fe 3d states which disperse only weakly in 
the /c^— direction.^ ^ Similarly to pnictide^^ the role of the 
dxy orbit can be approximated by a next nearest neighbor 
hybridization between d^z, dzy orbitals, and we consider a 
two-dimensional square lattice with two degenerate d^z, dzy 
orbitals per site. Based on these observations, our model 
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Hamiltonian is 
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k.cr 



where a is the spin index, r values are Pauli matrices and 

^ka ~ i'^xza ' '^yza] ^hc two Component field operator We 
consider up to third nearest-neighbor hopping for the present 
case as depicted in Fig. 1(b) which gives 
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Here c/sia — cos / sin (ifcao) for a = x,y, ji^ = fi 
and /i~ — 0. The Hamiltonian in Eq. 1 can be diago- 
nalized straightforwardly obtaining the eigenstates E^^ — 

fk ± \/(Ck)^ + (^k^)^-[Ref. 28] The TB parameters are ob- 
tained after fitting the dispersion to the LDA calculations,'^^ 

[Fig. 1(c)]. (ij., ty, txx: tyyi tx: ^y; ^xy) — ( 0.12, 

0.108, - 0.12, - 0.12, 0.048, 0.108, 0.06)eV and /i 
— 0.36eV. To obtain the experimental FS, we use a rigid 
band shift to /i = — 0.56eV and get two concentric squar- 
ish electron pockets centering (tt, 0), in good agreement with 
ARPES^ [see Fig. 1(d)]. 

As shown below, for such an FS, nesting along Q = (tt, tt) 
is most pronounced leading to a (ij.2_j,2— pairing according 
to conventional views. We include an SC gap with the same 
amplitude on each band to obtain the SC quasiparticle disper- 



sion as E^ — 



(-^ok)^ + where Ak — Ao(cosfc^a — 
cos kya) /2. Aq is taken from ARPES to be S.SmeVj^for both 
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FIG. 1. (Color online) (a) Unit cell of a KFe2Se2-layered com- 
pound, (b) x-y plane of the unfolded zone with Fe 2>dxz and idyz or- 
bitals. (c) The two resulting bands are plotted along high-symmetry 
momentum cuts and compared with the LDA one [dotted (black) 
lines]).' We mainly concentrate on fitting the low-energy part of 
the band to get the FS shape and their orbital characters in agreement 
with ARPES. (d) Computed electron pocket FSs overlain ARPES 
data^. (e) The normal state (NS) DOS is compared with its SC coun- 
terpart. An artificially large value of gap=100meV is chosen here to 
highlight the nodeless and isotropic nature of the gap opening at Ep, 
despite d— wave pairing. 



eigenstates throughout this paper, except in Fig. 1(e) where an 
artificially large value of Aq is chosen to explicate the behav- 
ior of unconventional gap opening on the DOS. Despise the 
presence of nodes in the underlying gap symmetry, the disper- 
sion at all momenta and the resulting DOS exhibit nodeless 
behavior in Fig. 1(e) as the FS pieces in these compounds 
are small in size and reside centering at the four-fold sym- 
metric momenta of the gap maxima, similarly to pnictide but 
unlike in cuprate. The nodeless and isotropic gap is observed 
in ARPES^Eand specific heat measurements^. 
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FIG. 2. (Color online) (a) Bare static xoil,'^ = 0), and (b,c) its 
RPA part, x'il)' for representative values of U. (d) x'il) for 
different values of U is drawn along the high- symmetry lines. The 
black and green arrows separate the peaks in x'il) coming from the 
inner and outer FS pockets (see text), (e) Schematic top view of the 
FeSe layer with the Fe spins in the stripe AFM order expected in our 
calculation are shown by the arrows. 



III. STATIC SUSCEPTIBILITY AND STRIPE GROUND 
STATE 

We now proceed to the spin-susceptibility, x, calculation of 
our model Hamiltonian. The calculation of the BCS x for a 
multiband system including matrix-element effects of all the 
orbitals is standard and can be found for example in Refs.l27}- 
l30l We consider the electron-electron correlation on the same 
Fe atom within the RPA formalism.^^ The terms in the RPA 
interaction vertex Us that are included in the present calcula- 
tion are the intraorbital interaction U, an interorbital interac- 
tion U' = U ~ J/2, the Hund's rule coupling J = [//4, and 
the pair hopping strength J' = jW^ 

Due to the specific nature of the FS shown in Fig. 1(d), 
the interpocket nesting is the dominant and incommensurate 
centering at q — (7r,7r), see static bare x'oil) in Fig- 2(a). 
An additional nesting occurs centered at g = due to the 
intrapocket nesting of the semi-squarish portion of the FSs. 
While the shape of XaiQ) is governed by the shape of the FS, 
the corresponding intensity is related to the Fermi velocity of 
the corresponding 'hot-spots'.^' Each of the x'q{q) pattern is 
split into two due to the splitting of the electron pockets into 
two concentric ones within our model Hamiltonian. 

Within the RPA, at different critical values of U, one or the 



other piece of x' (?) will dominate, leading to a SDW order 
phase, see Fig. 2(b)-2(e). To clarify that, we study the evolu- 
tion of the RPA susceptibility at different values of U, keep- 
ing the other interactions same with respect to the value of U 
(following the expressions given above). For a broad range 
of U studied in this case, the incommensurate peak around 
q — (tt, tt) is the winner, suggesting that a stripe like the 
SDW order will be present in this class of materials. These 
observations are consistent with the LDA calculations^ ^ as 



well as with number of experiments .Q^QIl In the small range 
of U, the inner electron pocket gains more intensity (as in the 
case of indicated by black arrows), while the outer one 
is comparatively small (green arrows). Interestingly, when 
U > 2.2eV, the situation is dramatically reversed and the 
outer nesting develops remarkable instability properties both 
around q — (tt, tt) and q ~ 0, although the former is always 
dominant in the range of U studied here. With careful obser- 
vation, one can find that the strong intensity within RPA is 
shifted toward q = (tt, 0.57r)/(0.57r, tt). This specific incom- 
mensurate vector is responsible for the possible 45° rotation 
of the dynamical spin resonance spectra studied below. 

The small q = (0, 0.467r) /(0.467r, 0) nesting which we find 
to be considerably large for our choice of parameters is close 
to the q = ferromagnetic (FM) instability or may induce 
some form of charge-density wave (CDW) with finite q mod- 
ulation. Note that the coexistence of a CDW and SDW due 
to multiple nestings is a n em ergent phenomenon predicted 
in pnictidJSl and cuprateP^SI The presence of vacancy or- 
der in these systems can p romote the coexistence of sev- 
eral competeting order.ESSl Scanning tunneling microscopy 
(STM) may be able to detect it. 

Unlike the (tt, 0)-CDW order in cuprates or the (tt, 0)- 
SDW order in pnictide, here the q — (tt, 0.57r)/(0.57r, tt)- 
stripe SDW order will facilitate d-wave pairing. This is due 
to the fact that the q — (tt, 0.57r)/(0.57r, 7r)-nesting comes 
from the inter-electron pockets at which the d-wave pairing 
possesses a change in sign of the gap, which is the criterion 
for spin resonance to occur as calculated below. But q ~ 
(0, 0.467r) / (0.467r, 0) instability will provide pair-breaking 
contributions to the d— wave symmetry. 



IV. DYNAMICAL SUSCEPTIBILITY, d-WAVE PAIRING 
AND SPIN RESONANCE 



Based on these findings, we proceed to study the evolution 
of the spin resonance spectra in the SC state, i.e., 1^) 
(RPA-BCS), which is directly measured by an inelastic neu- 
tron scattering (INS) experiment. In the SC state, x'o (Q; 
is strongly enhanced for ut < 2A due to turning on of the 
pair-scattering terms. This reduces the critical value of the 
interaction U to satisfy the RPA denominator to be positive 
at C/xg(q, < 1. We show the resonance spectrum for 
U = 1.5eV in Fig. 3(a) and that for U = 2.14eV in Fig. 3(b), 
which obtain a commensurate and an incommensurate reso- 
nance, respectively, in addition to their characteristic disper- 
sion. 

The dispersion of x"(q. which has a one-to-one corre- 
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FIG. 3. (Color online) Computed BCS-RPA x"{<l,^) for (a) U — 1.5eV and (b) U — 2.14eV. The constant energy spectra at two represen- 
tative energies are shown in the lower panels of the corresponding spectra. The right-hand column shows the neutron data on FeTei-j^Sej,, 
demonstrating the observations of (c) an upward dispersion'', in (d) an 'hour-glass' or 'camel-back' pattern'^ and (i,j) 45° rotation of the 
spectra^ in this compound. 



spondenc e with the dispersion of XoIQi comes from three 
condition J^mmnilll in addition to the matrix-element effects: 

(1) a..e.(q) = Ek-„k^; lAtc^l + |A^„-|, (2) sgn(Ak^) ^ 
sgn(Aj.„'), and (3) q = + where v, v' are the band 
indices. Condition (1) comes from the energy conservation 
of inelastic scattering of Bogoliubov quasiparticles, while 

(2) constitutes the coherence factors of BCS susceptibility 
to possess a non-zero value. As both FS pockets have the 
same gap symmetries and amplitudes, and particularly for 
FeSe systems, where all the FS pockets have isotropic gaps 
at all ki? values, conditions (1) and (2) are satisfied equally 
at all kp values for inter electron pocket scatterings with 
(i2;2_j,2— symmetry. No intensity modulation in xo('?) is gov- 
erned by conditions (1) and (2) and the only source of in- 
tensity variation is the inteiplay between the orbital matrix- 
element effect and condition (3), the latter is simply related 
to the number of degenerate ki? values. This means that the 
upward dispersion of the resonance behavior in FeSe can not 
be understood by tracking how the gap varies on the FS as 
happenes in cuprates,'^but only from the FS nesting and the 
matrix-element in the tensor form of and, also, on the in- 
teraction vertex Ug 



Furthermore, unlike in cuprates, where the 'hour-glass' dis- 
persion extends to w = due to the presence of nodal quasi- 
particles on the FS, in the nodeless FS of FeSe-compounds 
the resonance behavior spans only within A to 2A. For 
U = 1.5eV, the maximum intensity lies at the bottom of the 
spectrum at the commensurate vector at cj^es = 13.5meV 
with Wres/2A = 0.79. Again, U = 2.14eV, the maxi- 
mum intensity shifts to the top of the spectra at an incom- 
mensurate vector q ~ 0.78(7r, tt) at w^es — 12.4meV with 
a;res/2A = 0.73, slightly larger than its average value of 
0.64 found in all other superconductors.^^ The upward disper- 
sion of the magnetic spectra is qualitatively analogous to the 
one observed by Neutron in FeTeo.6Seo.4 superconductorsSl 
shown in Fig. 3(c), while pnictides show more commensurate 
resonance within the experimental resolution.^" By including 
the weak-intensity of the dispersion below the (7r,7r)- reso- 
nance, one can draw an 'hour-glass' or 'camel-back' feature, 
also seen in chalcogenideJ22lin Fig. 3(d). 

The constant energy profile of the resonance spectra is more 
interesting and bears important physical insights. Atw > ojres 
for both values of U, the resonance profile is squarish with 
maximum intensity lying along the diagonal and finite inten- 
sity along the (100) direction, see Figs. 3(f) and 3(h). For 
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FIG. 4. (Color online) Correspondence between (a) the 1 Fe unit cell and (b) 2 Fe unit cell in the FS topology and SC gap symmetry. The black 
to white color scale denotes the negative to positive value of the gap. (c) Transformation of the d^^y-wave gap and the FS are plotted in the same 
2 Fe unit cell as in (b) but with one Fe vacancy. (d,e) Real space view of the gap-symmetry corresponding to momentum space ones is given 
in (b) and (c). The dj.2_j^2-wave in the 1 Fe unit cell (magenta box) automatically takes the form of a d^y ~ sin {kxa/2) sin {kya/2)-wave 
in the 2 Fe unit cell (blue box), (e) Vacancy order obtained from the experimental data in Ref.|37| Yellow dashed lines in each plot denote the 
nodal lines for the corresponding pairing symmetry. 



U — 1.5meV, the intensity remains concentrated at (7r,7r) 
even when lj < ujres- Interestingly, for U — 2.14eV at 
Lo < LOres the peak along the diagonal sharply looses inten- 
sity while the intensity peak along the bond direction remains 
almost the same. As a result, the resonance profile is rotated 
by 45° as one goes down below the resonance. The 45° rota- 
tion of the resonance spectra is observed earlier in hole doped 
cuprates'^ and can also be seen in chalcogenides between two 
different dopings in the SC state^^ in Fig. 3(i) and 3(j). 



V. ONE FE UNIT CELL TO TWO FE UNIT CELL 
CONVERSION 

It has been argued that when FS pockets are translated 
from a 1 Fe unfolded zone to a 2 Fe folded zone, it switches 
location from k — (tt, 0)/(0, tt) to fc = (ti'iTt). Thus in 
a (ij.2_j^2— symmetry, the SC gap on the FS is assumed to 
change from nodeless to nodal behavior.^'^ We show that this 



apparent inconsistency is merely a result of misinterpretation 
coming from not performing the same unitary transformation 
in the gap symmetry as in the FS. To perform this transforma- 
tion, we denote a quantity with a tilde in the 2 Fe unit cell. 

The unitary transformation consists of 45° rotation of the 
lattice with a = ^/2a which gives k^/y — [kx ± ky)/2. 
In doing so, we find that the pairing symmetry (ij,2_j,2 ~ 
cos {kxo) — cos {kyo) in the 1 Fe unit cell transforms to an 
extended d^y ~ sin [k^a / 2) sm {kya / 2) , see Figs. 4(a) and 
4(b). This transformation can be more clearly understood 
from the real space view of the gap symmetry in the lattice 
in Fig. 4(d). The 1 Fe unit cell and 2 Fe unit cell are high- 
lighted by magenta and blue boxes, respectively. It is clear 
that the nodal line passes through the diagonal direction in the 
former case, which automatically becomes aligned along the 
zone boundary in the latter unit cell. The rotation of the SC 
gap phases follows similarly. We note that such unitary trans- 
formation preserves the nodeless and isotropic properties of 
the SC gap in any unit cell. 
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FIG. 5. (Color online) (a) FSs in the 2 Fe unit cell BZ are shown 
with various nesting channels (arrows), (c) Bare static susceptibility, 
Xo, is plotted in a 2 Fe unit cell, (c) Corresponding RPA instability at 
U — 1.7eV. The unitary transformation leads to both nesting vectors 
of the 1 Fe unit cell, shown in Fig. 2(a), coinciding around q ~ 
{n/2, n/2). (e) The dynamical spin suscebtibility, x" is calculated 
at U_— l.SeV as shown along the zone boundary [green arrow in (b)] 
for dxy pairing in the 2 Fe unit cell [dx2 _y2 in the 1 Fe unit cell], (e) 
Computed spin-susceptibility is plotted for nodal ci^2_j,2 -pairing in 
2 Fe unit cell [dxy in 2 Fe unit cell]. Here U is chosen to be 1.5 eV. 



An important consequence of the dxy-wave SC gap is that it 
breaks the crystal symmetry. As can be seen from the momen- 
tum space view of the pairing phases in Fig. 4(b), for example, 
as we translate from k = (— tt/o, —tt/q) to k ~ (tt/o, — 7r/a) 
by the reciprocal vector 27r/a, the phase of the SC gap is re- 
versed. One can recover the translational symmetry of the 
crystal by moving to fc = (27r/a, — tt/o) by a reciprocal vec- 
tor of 47r/a. In the real space approach, this would correspond 
to a case when one Fe atom from the center of the 2 Fe unit cell 
is removed; see Fig. 4(e). Remarkably, such a crystal symme- 
try can be achieved when the vacancy order and the magnetic 
order conspire with each other to favor the d^y-wave SC gap. 
This vacanc^arrangement was found indeed to occur in the 
real system.^"* In this context, we note that superconductivity 
occurs in this class of materials when a suitable amount of va- 
cancy is achieved which forms order and the cr ystal symmetry 
hops from the I4/mmm to the 1 4 /mm group .^3221 
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FIG. 6. (Color online) (a) Temperature dependence of the spin exci- 
tation spectrum at g — {n,n) for U = 1.5eV is compared with the 
similar observations by inelastic neutron scattering measurements (b) 
in pnictide^" and (c) chalcogeniges^*. Inset: Computed value of the 
intensity of the resonance peak is compared with BCS A(r). 



A. Rotation of spin resonance in 2Fe unit cell 

The same unitary transformation leads to both the nesting 
vectors of 1 Fe unit cell shown in Fig. 2(a) to coincide with 
each other around q = {tt/2,tt/2) in the 2 Fe unit cell as 
shown by blue arrow Figs. 5(a)-5(b). Subsequently the nesting 
along the (200)-direction shown by green arrow serves as the 
'hot-spot' to the d^y, where the SC gap changes signs. This 
leads to resonance spectra around {2tt ± 6, 0) /(O, 2tt ± S) as 
shown in Fig. 5(d). 

Note that the smaller intrapocket 'hot-spot' vector (blue ar- 
row) will give a sign change of the SC gap on the same FS 
pocket and thus a nodal SC gap [dx2_y2 in the 2 Fe/unit cell 
and dxy in the 1 Fe unit cell] will occur The resulting INS 
spectrum in a nodal SC state is confined near q ~ and does 
not show any true resonance peak but a spin-excitation dis- 
persion as shown in Fig. 5(e). In this case the associated weak 
magnetic excitation strength may not be sufficient to produce 
large values of Tc in these materials. 



VI. TEMPERATURE DEPENDENCE OF THE 
SPIN-RESONANCE AT Q 

In pnictide and chalcogenide superconductors the temper- 
ature evolution of the spin resonance follows the BCS mean- 
field behavior of the SC gap A(T)!^2|241. We predict that sim- 
ilar phenomenon also occurs in FeSe-superconductors. To 
prove that, we calculate the A(r) for the two band sys- 
tem solving the standard BCS gap equation with phenomeno- 
logical pairing strength parameters Vi,2= 52, 46 meV (ne- 
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FIG. 7. (Color online) Calculations with the five-band model in the 1 
Fe unit cellP^(a) Calculated band-structure is drawn along the high- 
symmetry lines for a chemical potential ^ = 7.45eV. (b) The cor- 
responding electron-pocket FS is plotted in the first quadrant of the 
Brillouin zone, (c)-(d) The bare and the RPA component of the static 
spin susceptibility, respectively, (e) Computed spin excitation spec- 
trum is shown along the diagonal direction in a d^2_y2— wave pair- 
ing symmetry. These results should be compared with the results 
obtained with the two band model in the 1 Fe unit cell in Figs. 1-3. 
The value of U for all RPA calculations here is set to be 1.5 eV. 



glecting interband pairing). The parameters are adjusted to 
obtain an experimental gap of Ao=8.5meV with BCS ra- 
tion 2A/A:bTc=6.7, in close agreement with the experimental 
value of 7.^ A(r) is shown in Fig 6(a) (cyan line) on top of the 
calculated RPA-BCS spin resonance spectrum at g = (tt, tt) 
with U = 1.5eV. The intensity at the resonance is plotted in 
inset to Fig. 6(a) in red line. Both the resonance energy and 
the intensity shows a remarkable one-to-one correspondence 
to the value of A(T). The results agree well with those for 
pnictide and chalcogenide shown in Figs. 6(b) and 6(c) re- 
spectively. The T— dependence of the resonance energy and 
its intensity follows the BCS form of A(r), suggesting that 
the itinerant RPA-BCS Fermi-liquid theory will be valid in 
FeSe-based superconductors in particular and in all iron based 
superconductors studied to date in general. 



but are essentially governed by the FS topology. The TB 
model for 5 bands is taken from Ref.i40i which includes the 5 
d orbitals of the Fe atoms in the 1 Fe unit cell. The band struc- 
ture is given in Fig. 7(a) and the corresponding FS is shown 
in Fig. 7(b). The topology of the FS is qualitatively captured 
here, although the five band model does not incorporate the 
nearly degenerate FSs as calculated in the LDA. Nevertheless, 
this minor technical drawback does not alter our conclusions. 

The magnetic structure including the two nesting vectors 
at q = (7r,0.57r)/(0.57r,7r) and g = (0, 0.467r)/(0.467r, 0) 
is well reproduced in the five band calculation as shown in 
Figs. 7(b) and 7(c). Furthermore, the magnetic resonance 
spectra is remarkably close to the one obtained in Fig. 3. The 
subtle differences between the two band and the five band 
calculation come from the inconsistency in the FS topology, 
not in the overall band dispersion. The remarkable consis- 
tency between the two band and the five band model suggests 
that the magnetic ground state as well as the superconducting 
properties can be studied properly once all the FS pieces are 
modeled accurately. 



VIII. CONCLUSION 

In summary, we have calculated the static and dynamical 
spin susceptibility for an effective two orbital TB model with 
d-superconductivity. We find that the leading instability in x' 
is evolved around an incommensurate vector q = (tTjO.Stt) 
at some critical value of U, pinning stripe like SDW or AFM 
order, in agreement with a number of experiments." Un- 
like in other iron based superconductors, here the SDW order 
has the same vector for the d-pairing symmetry to obtain 
the opposite sign at the 'hot-spot'. Our choice of parame- 
ters also predict a comparatively we ak C DW modulation with 
q = (0, 0.467r), similar to cuprates,-'' and pnictides.-'^ The 
d— wave gap yields a spin resonance behavior with an 'hour- 
glass' or upward dispersion and the 45° rotation of the reso- 
nance profile, in close agreement with observations in chalco- 
genides. The present results also demonstrate that the weak 
to moderate pair coupling theory of the spin fluctuation mech- 
anism of superconductivity is also appropriate in FeSe-based 
superconductors in particular and in all iron based supercon- 
ductors studied to date in general. A similar intermediate cou- 
pling theory has b een su ccessful in describing many salient 
features in cuprates.l^^J^ 

Recently, we became aware of the works by Wang et al^ 
and Maier et al^^ which also predicted the d— wave gap in 
FeSe based systems. 



VII. FIVE BAND MODEL 

We repeat the calculation presented in Figs. 1-3 within a 
five-band model and find that the magnetic structure and mag- 
netic resonance do not depend on the overall band structure 
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